Nonlinear thermoelectricity in point-contacts at pitch-off: a catastrophe aids cooling 
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We consider refrigeration and heat engine circuits based on the nonlinear thermoelectric response 
of point-contacts at pinch-off, allowing for electrostatic interaction effects. We show that a refriger- 
ator can cool to much lower temperatures than predicted by the thermoelectric figure of merit ZT 
(which is based on linear-response arguments) . The lowest achievable temperature has a discontinu- 
ity, called a fold catastrophe in mathematics, at a critical driving current I — I c . For I > I c one can 
in principle cool to absolute zero, when for I < I c the lowest temperature is about half the ambient 
temperature. Chargeless particles (typically phonons and photons) stop cooling at a temperature 
above absolute zero, and above a certain threshold also modify the discontinuity. More generally, 
we show that any system with a high figure of merit should have its nonlinear response analyzed; 
since its figure of merit gives little indication of its potential as a refrigerator or heat-engine. 



Introduction. Nanostructures often have thermo- 
electric responses, with electrical-currents causing heat- 
currents, and vice- versa PJH]. There have recently been 
a number of proposals for systems with large thermoelec- 
tric responses [5 12 which could have engineering appli- 
cations for efficient thermoelectric power-generation and 
refrigeration. In particular, it is hoped that they could 
efficiently cool electrons well below the temperature of 
standard cryostats [13 , which are increasingly inefficient 
at sub-Kelvin temperatures. 

Here we calculate the fully nonlinear thermoelectric 
response of a point-contact at pinch off. This system's 
linear (and nearly linear) thermoelectric response is well- 
studied experimentally |14l [T5] and theoretically fLU [TB1 - 
18 . However a refrigerator which cools to significantly 
below the temperature of its environment will typically 
be far outside the linear regime |131I19| ; since the temper- 
ature difference is not small compared with the average 
temperature. Here, we apply the nonlinear version of the 
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Figure 1: Heat-current J(Ti s i,7) through a point-contact 
when driven with a current /, for negligible phonon or pho- 
ton heating. Blue indicates cooling of the island in Fig. [2^,, 
while red indicates heating. The solid curve is the steady- 
state (J = 0), with the catastrophe at 7 C . The straight line is 
the maximum current, 7 ma x, corresponding to infinite bias. 



Landauer-Biittiker scattering theory [2QH2I] to thermo- 
electric heat-transport |25] through point-contacts. This 
shows that the dimensionless figure of merit, ZT, ceases 
to be an accurate measure of the thermoelectric response 
outside the linear regime. Electricity generation is worse 
than linear-response theory indicates, but refrigeration is 
better (achieving much lower temperatures than linear- 
response theory predicts). Indeed, the lowest temper- 
ature of the refrigerator is a discontinuous function of 
the electrical current. This discontinuity — a fold catas- 
trophe in mathematical language — occurs at a critical 
current I c , and is beneficial to refrigeration. For currents 
I < I c the refrigerator cannot cool below a finite temper- 
ature (of order half the ambient temperature for I — > J c ), 
while for I > I c it passes the catastrophe point and can 
in principle cool to absolute zero (see Fig. [T|). 

In practice, a thermoelectric device's quality is re- 
duced by heat flow carried by chargeless particles; usually 
phonons and photons. In the nonlinear regime, we show 
the weak phonon or photon effects do not significantly 
affect the catastrophe, I c , but do stop the cooling at a 
temperature above absolute zero. However at a critical 
value of the thermal transport due to chargeless particles, 
the catastrophe becomes a discontinuity in the derivative 
of the curve lowest temperature versus I. 

Nearly-linear analysis and its breakdown. The 
usual "nearly-linear" analysis [T] takes linear response 
theory plus a Joule heating term, and enables one to 
quantifies devices in terms of their dimensionless figure 
of merit ZT = GS 2 T/{O cl + 8 ph ), where T is the de- 
vice temperature, S is its Seebeck coefficient, G and 
e i are electrical and thermal conductances of electrons, 
while Gph is the thermal conductivity of chargeless exci- 
tations, principally phonons and photons. This nearly- 
linear analysis predicts electric power generation (when 
the island is heated) with an efficiency 



= VZT+l-l A Tp 
V VZT+1 + 1 V T isi 
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where T; s i and To are the island and environment tem- 
peratures. Typically, ZT is taken at the temperature ~ 
\ (To + Tsi). Carnot efficiency corresponds to ZT — > oo. 
For refrigeration, it predicts that the lowest achievable 
temperature, T m - mi is given by 



Tnin/To — 1 



\ZT. 



(2) 



The origin of Eq. Q is easily seen, and enlightens us 
about the breakdown of the nearly- linear analysis. One 
starts with the heat-current out of the island II 



J(T ish I) 



6+ (T - T isl ) 



(3) 



for Peltier coefficient difference n_ = (II2 —H%), sum of 
thermal conductances + = Q± + O2 + 6 P h and sum 
of electrical resistances R + = G7 1 + G7 1 . The first 
two terms are linear-response terms, while the last is the 
Joule heating. The steady-state curve, J = 0, shows that 
the lowest Ti s i is a quadratic function of I. The parabola's 
minimum is T m j n in Eq. Q — using the Onsager relation 
IT_ = S-T to note ZT = G + TP_/{Q+T). 

For our model of a point-contact at pinch-off (detailed 
below) , linear-response scattering theory [25 , 29 - 34J gives 
Gi = (e 2 /h) (1/2), III = -(fc B T /e)21n(2), and 9j = 
(k^To/h) (tt 2 /6 - 2[ln(2)] 2 ) [35,. This gives ZT ~ 1.4, 
which would imply that 



r/ = 0.22 (1 - T /T is 



T„ 



0.3 T . (4) 



However, this analysis fails when nonlinear effects are too 
strong to start from Eq. ([3| . Typically this occurs when 
the leading nonlinear Peltier term [T7], of the form II/ 2 , 
is larger than the Joule heating term ^R + I 2 . For the 

point-contact this is so, since II = (h/e 2 ) 3 In [2]. Then 
including this term in Eq. |3| changes the sign of the 
prefactor on I 2 ; which means we must look beyond I 2 to 
find the steady-state curve's minimum. Our fully nonlin- 
ear analysis below, shows that expanding about equilib- 
rium to any order |17| will not give this minimum. 

Nonlinear analysis. The thermoelectricity literature 
discusses J(Tj s i, I) — as in Eq. (J3J above — rather 
than J(X| s i, V) for voltage drop, V. This is because differ- 
ent thermoelectric devices are arranged in series electri- 
cally (see Fig.[2^,,b), so / is the same in all of them (unlike 
voltage drops). Thus it is easier to get response of a se- 
ries of elements from each element's J(T; s i,7) than from 
each element's J(Tj s i, V). For complicated non-linear re- 
sponses, the former is straight-forward while the latter is 
extremely difficult; thus we consider J(Tj s i,/). 

We take the island to be classical; i.e. big enough for 
particles entering it to thermalize to a Fermi distribution 
at temperature T; s i before escaping. See Refs. |26lf2"5] 
for cases where there is a quantum dot in place of the 
classical island. In our case, each point-contact can be 
treated by a separate Landauer-Biittiker scattering ma- 
trix analysis [25], see also Refs. |29If3"4"] . We generalize 
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Figure 2: Thermoelectric circuits made with point-contacts 
shown in (a,b); "e" ("h") means the point-contact is in a ma- 
terial whose charge-carriers are electrons (holes). One can 
minimise the heat-current carried by phonons and photons, 
Jph, by suspending the island |39j . (c) Motion of charges in 
the gates (arrows) caused by making Vl positive, which par- 
tially screens Vl at some distance from the point-contact, (d) 
Point-contact 1 tuned to pinch-off (E pc equals the island's 
chemical potential) by adjusting V g . 



these heat currents to the nonlinear regime |36| . includ- 
ing electrostatic (Hartree-like) interaction effects in a self- 
consistent and gauge-invariant manner; as Refs. |21If2"4"] 
did for charge-current, see also |37l [35]. To go beyond 
the voltage-squared contributions to transport (which 
Ref. [2T] treated in detail), we use a simple model of 
interactions, which is none the less gauge-invariant and 
self-consistent. The heat-current into lead i of a given 
nanostructure is 



Ji 



de 



£>-<M)-M{e-«w) m (5) 



where fj(e) = (l + exp [(e-qVj)/(kBTj)]) is the Fermi 
function, and q is the charge of the carriers; electrons with 
q = —e in point-contact 1 and holes with q = e in point- 
contact 2. The energy e and all voltages 14 are measured 
from the same external reference. The formula for the 
charge-current Ii into lead % is the same but with q in 
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place of (e — qVi) as the first factor in the integrand. The 
transmission function of a particle through the nanos- 
tructure from lead j to lead i is As 



Tr 



USij-S^de-qVk}) Sade-qVk}) 



t3 ({e - qV k }) 
where § 



the scattering matrix from lead j to lead i, and the 
trace is over all modes of those leads. Here Sy must be 
found self- consistently; it depends on the charge distribu- 
tion in the nanostructure, which in turn depends on Sjj. 
By writing as a function only of energy differences, 
{e — qVk}, we make the gauge-invariance explicit; i.e. it 
always satisfies [21] [(d/de) + J2k {d/ A(qV k ))]A l] = 0. 

Point-contact 1 is a two-lead nanostructure with elec- 
tron charge-carriers (q = —e). Having established the 
gauge-invariance, we are free to measure all energies e 
and all voltages V k (including V g ) from the chemical po- 
tential of the island (the point-contact's M lead). When 
Vl is non-zero, we assume that a proportion (1 — a) of 
this bias is screened by the electrostatic gates a long 
way from the narrowest-point of the point-contact, while 
the rest is screened self-consistently by the electron-gas 
(Fig. [2ji) close to the point contact. Then the screened 
point-contact induces a potential barrier of height, E pc 
(measured from the island's chemical potential), typi- 



cally obeying E pc - 
screening. Here E. 



E s = g scr (aqVh), where g sct is due to 
can be tuned at will, since it is eV E 
minus a geometry-dependent constant. Assuming a long 
enough point-contact that there is negligible tunnelling, 
one has j4lm(c — E pc > 0) = — 1 (perfect transmission) 
and Aim(c — -Epc < 0) = (no transmission) (3D]- As 
an example, the Supplementary Material gives a simple 
model of screening for which we derive g S cr( a '?^L) self- 
consistently. However, in what follows we allow the na- 
ture of screening (both a and the form of g SC r(aqVL)) to 
be completely arbitrary. Gauge-invariance is satisfied, as 
seen by noting that e — E pc is a function of {e — qV k }, 
since VL ig are biases relative to the M lead. 

To operate at pinch-off, we tune V g so that E pc = for 
any given Vl ■ If the electrostatic gates dominate screen- 
ing (a — ¥ 0), then E pc is independent of Vl, making it 
easy to tune to pinch-off. Otherwise the point-contact 
should be calibrated prior to use; finding the pinch-off 
point (the V g at which current starts to flow), as a func- 
tion of Vl- Then the charge and heat-currents from 
point-contact 1 into the metal island are 

/(2k, VL) = ^ [r M ln(2) - To In f 1 + e -^/k»T Y\ ^ 



Ji(2k,V L ) 



h 2 



2i^+ToLi 2 (- 
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for a dilogarithm function Li 2 (z) = J °° tdt/(e t /z — 1). 
From Eqs. ^f[ , we get 



Ji(2k,7) 



k 2 T 2 



12T 2 



Li 2 (l-exp p(T M ,J); 



,(8) 



where we define 3 = ft.[imax(2isi) — f\l (eksTo) and note 
that I < P nax (T is i) = eyfcBT is iln[2]//i. This function is 
given by the color plot in Fig. [l] For point-contact 2 
(in which the carriers are holes not electrons) we take 
-e «-> e, then J 2 (T ish I) = Ji(T ish I) since I 2 = -I. 
For J <C 1, we can use Li 2 (2:) = z + 0[z 2 ] to write 

Ji = (klT*/ti) [J - (7r 2 /12)(T isl /r ) 2 + 0[J 2 ]] , (9) 

so Ji as a quadratic in temperature and linear in current; 
the reverse of the nearly-linear theory in Eq. ([3]). This 
approximation captures the features of the exact result 
plotted in Fig.[T] except the top-left corner. This corner is 
the linear-response regime (small (T — T- ls \) and I), where 
one has Eq. Q with Li 2 (-1 + z) ~ -tt 2 /12 + ln[2]z. 

Refrigeration neglecting phonons or photons. 
Heat flow into the island is Jtotai °c J\ for the devices in 
Fig. (2^,,b; Jtotai = 2Ji for the thermocouple. The black 
curves in Fig. [T] are Jtotai = 0, giving the steady-state 
temperature (solid for stable steady-states and dashed 
for unstable ones). Solid curves give the temperature the 
island will be cooled to by a current /. Eq. (|9| tells us 
the steady-state has J as a quadratic function of T; s i; this 
approximation gives the catastrophe at eI c /(ekB,Ta) — 
3(ln[2]/7r) 2 ~ 0.14 with T isl /T = 61n[2]/7r 2 ~ 0.42, 
which is very close to the exact solution in Fig. [T| 

Refrigeration with phonons or photons. Assume 
that chargeless particles carry heat ballistically from hot 
to cold (a reasonable assumption for phonons in clean 
samples below afew Kelvins [41]), then J p h — alT^y—T^) 
where a — Aeo~sB for an island with surface area A and 
emissivity e, where ctsb is the Stefan-Boltzmann constant 
for the photons or phonons. Fig. [3^ shows that this moves 
the catastrophe to higher I and lower T; s i, until at a x 
/i(Tg/fce) 2 ~ 1/8 it is replaced by a discontinuity in the 
derivative of the steady-state curve. The effect remains 
qualitatively unchanged, if we take J p h = ch(Tq — T 2 sl ) as 
for some photonic transmission ]42_. 

The transition to a continuous response. A tran- 
sition from the response in Fig |] to a continous one 
(evolving towards the nearly-linear response in Eq. (|4l 



occurs upon reducing the e-dependence of Aij. This 
reduces the thermoelectric response compared to the 
normal electrical response. We do this by putting the 
point-contact in parallel with another system with e- 
independent transmission; e.g. a tunnel-barrier. Fig. [SJd 
shows the steady-state response of the point-contact in 
parallel with a load with conductance Gi oa d ~ (e 2 /h) x g 
for various g. For g > 1, the usual nearly-linear theory 
works. However for g — 1/4 we see the deviations (cf. 
solid and dashed curves) predicted by the rule of thumb 
below Eq. ^ , since n of the combined system is larger 
than i/iwhen g < (31n[2]-l)/2 ~ 0.53 03]. Atransition 
in the steady-state's behaviour occurs at g = g c ~ 1/200. 
For g < g c , the curve ceases to be single- valued, so T m j n 
becomes a discontinuous function of I. 
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Figure 3: (a) Steady-state refrigeration curves for different couplings to phonon or photons, «,. for a x hiTo/ks,) 2 = 
0, 1/32, 1/16, 1/8/, 1/4, 1/2, 1. (b) Steady-state refrigeration curves for a point-contact in parallel with a load with conductance 
Gioad = {e 2 /h)xg for g = 0,1/1000, 1/200, 1/40, 1/10, 1/4, 1. The solid curves are the exact results for J (T iBh I) = 0. The dashed 
parabolas are the nearly-linear approximation for g — 0, 1/4; it fails completely for g = 0, but improves with increasing g (for 
g = 1 there is no visible difference from the exact result), (c) Heat-engine efficiency curves for different strengths of coupling to 
ballistic phonon or photons. The different curves have coupling a chosen such that a x h(To/kB) 2 = 0, 1/256, 1/64, 1/16, 1/4, 1. 
The two dashed lines are the linear-response predictions for ah(To/kB) 2 = 0, 1. 



Heat-engine efficiency. For Tj s i > To, the circuit 
in Fig. [2^l provides electrical power P = IV to any load 
connected between L and R. Fig. [3J; shows the ratio of 
the optimized-power (load chosen to maximize P) to the 
heat-current; a nonlinear method for finding this is out- 
lined in the Supplementary Material. 

For Xisi 3> Tq we find that both the heat-flow carried by 
electrons and the electrical power go like (Tlsi/Xo) 2 , and 
the optimal efficiency saturates at 1 — yfl — 6(ln[2]/7r) 2 ~ 
15.9% in the absence of phonons or photons. We have 
no simple argument why this curve is slightly peaked at 
(Tfci - T )/T is i ~ 0.88. For finite coupling to ballistic 
phonons or photons, the curves decays to zero at large 
Tj s i, because J ph goes like (Ti sl /T ) 4 . If instead J ph — > 
a(Ti s \/To) 2 the optimal efficiency saturates at 1 



1 - 61n 2 [2]/(7r 2 + 6a), less than Eq. (l! predicts. 
Concluding remark. The principle experimental 
challenges to observing the discontinuity in Fig. [l] are 
minimizing phonon and photon effects, and keeping the 
point-contact at pinch-off. However, irrespective of 
whether this is easily achievable, our results give a im- 
portant message for all works on systems [SHU] show- 
ing large figures of merit, ZT. Having ZT > 1 is a 
strong hint that the nonlinear Peltier coefficient n may 
be larger than ^R. If so, optimal response will occur 
deep in the non-linear regime, where ZT ceases to give 
insight into the device's true potential as a refrigerator 
or heat-engine. This is already the case for the point- 
contact whose ZT is only 1.4; so we can certainly expect 
the same for systems with ZT ^> 1. In some cases (such 
as the point-contacts studied here) the device may be 
much better than its ZT would imply. 
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SUPPLEMENTARY MATERIAL 

Self-consistent solution for linear screening. We 

model the point-contact as a one-dimensional scattering 
problem (along the x-axis), with the potential sketched 
in Fig.[2jl. Close to the point-contact, this takes the form 
qV(x) — E s — kx 2 + qV SCI (x — x pc ) with energy measured 
from the chemical potential of the island. The trans- 
verse confinement generates the (Eg — Kx 2 )-term., where 
E g can be tuned, since it equals eV g minus a geometry- 
dependent constant. The gV^cr-term is the screening due 
to the electron gas, which we take to be of the form 



V acr (x) = < oVlQs 



ciVl for x < —l SCI 

- x)/(2l SCI ) for |af| < Z scr 
for x > l SCI 



with x pc self-consistently determined peak of qV(x). 
A little algebra gives a; pc = — aqV^ / (AkI sct ) , thus 
the energy at the peak is E pc = qV(x pc ) = E g + 
\aqV-L (l — a 9^ / L/(8K/g CI .)) . Finally we note that both a 
and Z scr depend on the scattering matrix of the junc- 
tion, which in turn depends on E pc . To solve this prob- 
lem self-consistently, we assume we are in the regime 



where e 



pc 



E„ 



j pc — Eg is small enough to approximate 
a = a (l + b a e pc ) and Z scr = l scl0 (l + he pc ). If necessary, 
«0! 4cr0i b a , bi can be found by simulating Poisson's equa- 
tion; typically e pc is small for small a. Then E pc is equal 
to a linear function of itself; re-arranging this gives 



E pc = Eg 



a Q qV L /2 - C{qV L ) 



1 - a Q b a qV L + 2C{qV h ) [b a - h 



where we define C(qV\ J ) — (a qVh / l SC io) /(16k). As 
mentioned earlier, we assume that tunnelling at ener- 
gies e < E pc is negligible, so Aua{c - E pc > 0) = -1 
and A-^m^ — E pc < 0) = 0. To see that this respects 
gauge-independence, we recall that e,-Epc,VL ig are all 
measured relative to the island's potential, and replace 
them by quantities measured from a fixed external ref- 
erence, so the island is at Vm- For clarity, here (un- 
like in the paragraph containing Eq. ([5])) it is necessary 
to use a tilde to explicitly indicate quantities measured 
from the external reference. We make the replacement 
qVh — (e — qVu) — (? — qVh). From this, we see that 
•Alm(e — Epc) is only a function of the set of differences 
{e^— qVk}, and so respects gauge-invariance. 

Voltage and power dependence of cooling. Ear- 
lier we explained that it is best to consider the nonlinear 
response of thermoelectric devices as a function of cur- 
rent. However, for completeness we give the response of 
a single point-contact as a function of voltage or power. 



In Fig. SI i the dashed-lines are contours of V(I, Ti 8 \) 
and solid-curves are contours of P(J, T isl ) = IV {I, T- m {); 
both V and P increase from top to bottom (with infinite 
voltage and power on the boundary of the red-triangle). 
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Figure SI: (a) Reproduction of Fig. [I] with lines of constant 
voltage (dashed) and constant power (solid) superimposed, 
(b) heat-current as a function of V and Ti s i, with the steady- 
state, J = 0, marked by the black-curve. The color-scales are 
as in Fig.[l] but with white used for all hJ/(k B T ) 2 < -0.25, 
and the darkest color used for all hj j '(ksTo) 2 > 0.18. 



These curves show that the steady-state is a continuous 
function of V, but a discontinuous function of power; the 
former is clearly seen in the plot of Eq. in Fig. |ST 



However one cannot extract the voltage-dependence of 
cooling for a circuit from this J(Xl s i, V), because the rel- 
ative voltage drop across each thermoelectric element de- 
pends nonlinearly on T[ s \. In general, the only way to get 
that response is to use Jtotai = Y]j Jj{Ti s i,I) where we 
sum j over each thermoelectric element. One can then 
use V(Ti s \,I) for each element, to find the total voltage 
drop that generates / at a given T^. 

Heat-engine efficiency in the non-linear regime. 
To calculate the maximum electrical power P(Ti s \, I) that 
a heat engine can extract from a heat flow J(T s i, I), one 
assumes a Ohmic load — so U(T s i,I) = I/G\ oa d — is 
connected across its terminals, and adjust Gi oa d to op- 



timize the ratio of the power at the load P(Ti s \,I) — 
I V(Ti s \, I) to the heat flow J(Ti s \, /). This corresponds to 
finding the / = I opt which maximizes P(T; s i, /)/J(T s i, I). 
The saddlepoint of P/J is given by P'J = PJ' 
where primed indicates (d/dl), If we have U(Tj S i,I) 
and J(Ti S \,I) we can solve this to find the optimal 
value of the charge-current I op t- The optimal efficiency 
is r\ — P{I op x)j J[I op t)\ although one also has 77 = 
P' (lopt) / J' (lopt) 1 which is often easier to evaluate. 

As a warm-up, we consider the usual linear problem, 
with V(T ish I) = S(T isl - T ) - G- 1 / and J(T U ,I) = 
6 (7k - T ) + III, with n = Y/G and S = B/G. 
Using these equations to calculate the optimal effi- 
ciency in the manner described above, we find / op t = 
(0/n)[v^(3k)2k + l - 1] (T M - T ). Assuming ZT is 
approximately T-independent, we get Eq. (JT|). 

Now we use the same method to get the efficiency in the 
nonlinear regime. Unfortunately, in general we cannot 
get an analytic solution of P'J = PJ' from Eqs. (6j8l, so 
we solve it numerically to find / opt for different Tj S i/To, 
and plot 77 against (T is i — T )/T is i in Fig. [3J;. However 
we can get an analytic result for large (Tj S i/To), using 
the fact (confirmed by the numerics) that in this limit 
-eUi(/ pt)/(fcBTo) > 1 with Vi < 0. Using ln[l + e"] -> 
H and Li2(— e M ) — > — |/it 2 for large /i, we get 
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where 

2 



we define t; s ] = Ti s \/To and k = 7T /6 — 
ln^[2] ~ 1.16. The heat-current from the hot source 
into the device is J(T; s i,7) = — Ji(T; s i,/), and P = 
-VtPisi, 1)1 (given that Vi < 0). In this case 
P'(T iah I opt )J{T iah I opt ) = P(T ish I opt )J'(T ish I opt ) is a 
quadratic equation for I pt, which we solve to find 



hl 



pt 



1 + \u[2]/k - 1 



Tisi 



ek B T ln[2] 

from which we get the optimal efficiency given in the text. 



